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Abstract 



Dirac equation is solved in the near horizon hmit geometry of an extreme 
Kerr black hole. We decouple equations first as usual, into an axial and 
angular part. The axial equation turns out to be independent of the mass 
and is solved exactly. The angular equation reduces, in the massless case, to a 
confluent Heun equation. In general for nonzero mass, the angular equation is 
expressible at best, as a set of coupled first order differential equations apt for 
numerical investigation. The axial potentials corresponding to the associated 
Schrodinger-type equations and their conserved currents are found. Finally, 
based on our solution, we verify in a similar way the absence of superradiance 
for Dirac particles in the near horizon, a result which is well-known within 
the context of general Kerr background. 



I. INTRODUCTION 



During the last three decades the study of spin-| particles on type-D spacetimes has 
attracted much interest and by now accumulated results are already available in the literature 
[1-3] (and references cited therein). The main reason for this is that all well-known black 
holes (BHs) are in this category and their better understanding involves detailed analysis of 
various physical fields in their vicinity. Dirac particles with (and without) mass constitutes 
one such potential candidate whose interaction/behavior around BHs may reveal information 
of much significance. Test Dirac equation in spacetimes other than BHs are also started to 
arouse interests for various reasons. Prom this token we cite the Robertson- Walker and 
Bertotti- Robinson (BR) spacetimes [4-6]. The latter in particular has already gained much 
recognition in connection with extremal BHs, higher dimensions and the brane world. 

In this paper we consider the Dirac equation in the near horizon geometry of an extreme 
Kerr BH. This constitutes the most important region of the outerworld prior to the horizon 
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of a Kerr BH enhanced with the extremahty condition. Extremal BHs are beheved to have 
connection with the ground states of quantum gravity. This alone justifies, in spite of the 
absence of backreaction effects, the importance of spin-i particles on such backgrounds. 
The throat geometry is a completely regular vacuum solution with an enhanced symmetry 
group SL{2, R) x U{1). In many aspects, this solution shares the common features with the 
AdS2 X 5*2 geometry arising in the near horizon limit of extreme Reissner- Nordstrom BH. 
The behavior of massless scalar fields in the extreme Kerr throat has been considered, and 
it is found that certain modes with large azimuthal quantum number exhibits superradiance 
[7]. This implies that the geodesies near the horizon can escape to infinity carrying energy- 
momentum more than the amount that infalls. Our solution enables us to investigate a 
similar phenomenon with Dirac fields which turns out to be negative as far as superradiance 
is conserved. This result is in accord with the treatment of Dirac equation in the general Kerr 
background in the absence of an exact solution [1]. The analysis of fermions in Kerr- Newman 
background and absence of superradiance was shown first by Lee [3]. Our aim is to revisit 
this item - not only by separating equations and deducing results on general grounds - but 
rather obtaining exact solutions and employing them. The advantage of confining ourselves 
to the near horizon alone bears its fruits by allowing solutions expressible in terms of known 
polynomials. In the general background of Kerr family of BHs even this much remains a 
problem beyond technical reach. Meanwhile, it is important to note that in contrast to the 
general Kerr metric, the throat metric is not asymptotically fiat. This difference naturally 
shows itself in the potentials, too. In other words, the extreme Kerr throat metric represents 
such a local region that the geometry of interest is different than the geometry of the general 
Kerr metric. The steeply rising potential prevents any particle (or field) flow to inflnity to 
make superradiance. Hence, we can mentally say that the behavior of particles in two 
geometries must be considered separately. 

In order to separate equations we employ the well-known method due to Chandrasekhar 
so that we prefer to label the set of equations as Chandrasekhar- Dirac (CD) equations. 
We separate the y (axial) and 9 (angular) dependence in such a way that the resulting 
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axial equation remains independent of mass. This leads us to an exact solution irrespective 
of mass. The angular equation on the other hand depends strictly on the mass. For the 
massless case (which we refer to as neutrino equation) the angular equation reduces to a 
confluent Heun equation [8]. When the mass is nonzero, however, we can not identify our 
equations but instead we express them as a set of linear equations suitable for numerical 
analysis. 

Organization of the paper is as follows: in section II we review the near horizon geometry 
of an extreme Kerr BH and separation of variables of the CD equation. Solution of the axial 
equation follows in section III. The massless and massive cases are discussed in sections IV 
and V, respectively. The reduction of our equations into one-dimensional Schrodinger-type 
equations with their conserved currents and superradiance are all included in section VI. 
The paper ends with a conclusion in section VII. 

II. EXTREME KERR THROAT GEOMETRY AND SEPARATION OF DIRAC 

EQUATION ON IT 

The extreme Kerr metric in the Boyer-Lindquist coordinates is given by 



e^^ (d^ + u;dt)'-p'(^ + de' 




(1) 



where 



2v 



Ap2 




(2) 



Ap2 
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,2 ^ r'^ + M'^cos^e 



In the extreme case, both the total mass M and the rotation parameter a become identical 
so that the angular momentum J = and the extremal horizon corresponds to f = M. 
The area of the horizon is yl = Stt J. 

To describe the near horizon (or throat) limit of the extreme Kerr metric, due originally 
to Bardeen and Horowitz [7], one can set 

r = M + Ar 



(3) 



2AM 



and takes the limit A — > 0. In these new coordinates, the throat metric is obtained as 



^dr' - ride' 



Tg sin^ 9 



+ ^dt] 



(4) 



where 



F = 



1 + cos^ 



rl = 2M2 



(5) 



We set further, for simplicity, Tq = 1 . This throat spacetime has no longer asymptotic 
flatness. 

Finally, passing to more general coordinates, 

1 



{i + r')-i 



cot t 



2tr 



(6) 
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(f) = if + ln 



cos t + y sin t 



1 + tr 



we can write the throat metric (4) as follows 



\ J 1 + 1/ 



{d^ + ydtf 



(7) 



The metric functions in (7) depend only on the variable Q thus as expected in the search 
of solution for the Dirac equation the angular equation forms the crux of the problem. The 
coordinates — oo < f < oo , — oo < y < oo cover the entire, singularity free spacctime. 
The Killing vector ^ is not timelike everywhere; it admits a region (for sin^ 9 > 0.536) in 
which it becomes spacelike. Therefore by a coordinate transformation this particular region 
is transformable to the spacetime of colliding plane waves [9]. Recently, it has also been 
shown that the metric (7) can be obtained as a solution to dilaton-axion gravity which is 
similar to the rotating BR spacetime [10]. 

The singularity free character can best be seen by checking the Weyl scalar \1'2 and the 
Kretschmann scalar: 

2 



*2 = 



(1 + cos2 9y 



3 cos 9 — 1 + i cos 9 (cos^ ^ - 3) 



192 sin^^ 



(l + cos^ 9)^ - l&cos^9 



(8) 



(9) 



(1 + cos2 9f 

When the backreaction of the spin-^ test particles on the background geometry is ne- 
glected, the Dirac field equation is given by the CD equations [1] on a fixed space-time 
(7). 

We choose a complex null tetrad {Z,??., m,m} that satisfies the orthogonality conditions, 
l.n = —m.m = 1. We note that, throughout the paper, a bar over a quantity denotes complex 
conjugation. Thus the covariant one-forms can be written as 



V2I = ^F(l + y^)dt 



}dy 



V2\ 



n 



F{l + y^)dt + 



^dy 



(10) 



/- iysin^ , rr^ ..^ ism.9 , 
V2m = ^ dt + VFd0 H ^d(p 

r- iysm9 , /r; ,^ ism.9 , 

^^F a/F 



and their corresponding directional derivatives are 



If{1 + y2) x/F 



yF(iT 



2^ VF 



v^5 = -4^a,-^a^ (11) 



V2-5 = -^do + '^d, 
VF smt^ 



One can determine the nonzero NP complex spin coefficients [11] in the above null-tetrad 

as 

sin 9 (cos 9 — i) 



TT = — r = 



(2F) 



£-7= ; ' (12) 
2J2F(l + y2) 



2 cot 9 — i sin ^ 

a = — p = 



2(2F)^ 

The CD equations in the NP formalism are then[l] 



{D + e- p)Fi+ (5 + n-a)F2^ ifXpGi 
{A + 1^- ^) F2 + {5 + P - t) F, = t^pG2 
{D + e-p)G2-{5 + W-a)Gi^ilipF2 (13) 
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{A + ji--f)Gi- [5 + (3-T)G2^ifipFi 

where //* = V^fJ,p is the mass of the Dirac particle. 

The form of the CD equations suggests that we introduce [1,6] 

i^2 = /2(2/)A3(^)e^(^*+"^'^) (14) 
^2 = 52(^)^4 (^)e^('^*+"^'^) 



where a is the frequency of the corresponding Compton wave of the Dirac particle and m 
is the azimuthal quantum number of the wave. Our convention is that a is always positive. 

Inserting for the appropriate spin coefficients (12) with the spinors (14) into the four 
coupled CD equations (13), we obtain 



/2 "/2A1 A, 



^^V ■ .92 A, r- {L+A,] 
= -IfX -r — vF - - — - 

h h Az A3 



(15) 

91 91^4 A4 



,,'li^^.(£AA 

92 92 A2 A2 

where the axial and the angular operators are 

1 



1 



and 



Z=^l + y^dy + ^^7^ [y - 2^cT] (16) 



_ mF 1 f cos^ 9 ism 9' 
L^d0 + + — + — — 

sm9 2F \ sm9 2 



^, _ mF 1 f cos^ 9 ism9^ 

^ de \ \ 

sin^ 2F \sin9 2 



mF 1 f cos^ 9 ism9\ 



"^~^^2Fl^^i^ 2~J 
respectively. One can easily see that L = £ and L"*" = £+. 

Further, choosing fi — g2, f2 — Qi, — ^2 and A3 — and introducing the separation 
constant as iA, where A is a real constant, we can separate Dirac equation (15) into axial 
and angular parts 

Zgi = -i\g2 (18) 
Zg2 = iXgi (19) 

and 

LA3 + ^^l*A2^/F = iAAi (20) 

L+Ai + iii*Ai\fF - zAAa (21) 
It is clear from equations (18,19) that gi = '§2- 
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III. SOLUTION OF THE AXIAL EQUATION 



If we decouple the axial equations (18,19) in equation (18) to get gi, we obtain 

z(Zg,) = X'g, (22) 

Similarly one can decouple the axial equations in equation (18) for g2. The explicit form 
of equation (22) can be obtained as 

(l + y') g'iiv) + '^yg'M + 3^ Q + + ^ - (1 + y') + iay^ g,{y) = (23) 

( Throughout the paper, a prime denotes the derivative with respect to its argument.) 
Thus the solutions of the decoupled equations for g^i, equation (23), and g2 (not given 
here) can be found in terms of the associated Legendre functions as follows 

9i{y) = ci-Pfi {ty) 



gM = C2Pf_i {-ly) (24) 

^ 2 

where 



+ (25) 

and ci, C2 are complex constants. 

Here, due to the physical necessities, we considered only the first kind of the associated 
Legendre functions. Although solutions (24) seem like complex solutions, it is possible to 
draw the real functions from the above associated Legendre functions. We may define 

A = m+- with m = l,2,3.... (26) 
2 



and 



a^j (Vl6n2 - 3 - l) (27) 
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so that 

13 — n with n = — m, — m + 1, — 1,1, m — l,m (28) 

In order to get the real functions for solutions (22), the required condition is: m — 
\n\ =even number. 

It is worth also drawing attention to the following remarks: 

i) In the case of A = 0, equations (18,19) reduce to simple first order differential equations 
which admit the following solutions 

^i(y)=C3(l + ?/^)"^e-*--^(^) 

^2(2/)=C4(l + i/2)"^e--*--^(^) (29) 
with C3, C4 complex constants. 

These two solutions can be interpreted as representing ingoing and outgoing waves. 
a) In the case of A = |, we obtain the following complex solutions from equations (18,19) 

..(y)-.fi^)Vcsfl^)"^ (30) 
where again Cj with j = 5, 6, 7, 8 are complex constants. 

IV. REDUCTION OF THE ANGULAR EQUATION TO HEUN EQUATION: THE 

MASSLESS CASE 

In this section, we shall show that the angular equations (20,21) for the neutrino particles 
can be decoupled to the confluent Heun equation. To the end that let us reconsider equations 
(20,21) in the exphcit form for //* = 0, 

A'^{e) + {K + G)A3{e) = iXAi{e) (31) 
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A[{e) + {K- G)Ai{0) = iXAs{9) (32) 

where 

G = ^ (34) 

By introducing the scahngs 

A^{9) = Hi{e)e-I^''-^^'^ (35) 

As{9)^Hs{9)e-I^''+^^'^ (36) 



one gets 



H[{9)^iXHs{9)e-I^^'^ (37) 



h;,{9) = iXHi{9)eI^^'^ (38) 
If we decouple equations (37,38) in equation (37) for Hi{9), we get 

H';{9) + 2GH[{9) + X^Hi{9) = (39) 
In similar fashion, we find, for H^{9), 

H'^{9) - 2GH'^{9) + X^H^{9) = (40) 
Introducing a new variable 9 — cos~^ (1 — 2z) , equations (39,40) turn out to be 

H';{z) + (-2m + + H[{z) - ^^^Mz) = (41) 

miz) + (-2m + ^ + ^) H!,{z) - ^^^H,{z) = (42) 

Let us recall the general confluent form of Heun equation [8] , 
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Drawing the similarities between equation (43) and equations (41,42), we observe the 
following correspondences: 

a) For equation (41), 

D^Q, X^, A = -2m, 5=^+m and C = ^- m (44) 

b) For equation (42), 

D^O, /i = A^ A = -2m, S = m- ^ and C = m+ ^ (45) 

The confluent Heun equation (43), with its accessory parameter h, has two regular sin- 
gular points at 2; = 0, 1 with exponents (0, 1 — B) and (0, 1 — C), respectively, as well as an 
irregular singularity at the infinity. In the vicinity of the point z = 0, its power series can 
be written as 

00 

H{D,A,B,C,h;z) ^Y^WjZ^ (46) 

j=o 

and the coefficient Wj satisfles a three-term recurrence relation [8] , 

Wo = 1, ^1 = ^ (47) 

(j + 1) (j + B) - A{j - 1 + D)W,., = [j{j-l-A + B + C)- h] W, 

It is also possible to obtain the power series solution in the vicinity of the point 2; = 1 by 
a linear transformation interchanging the regular singular points z — and z = 1. Namely, 
z ^ 1 — z. 

Expansion of solutions to the confluent Heun equation in terms of the hypergeometric 
and confluent hypergeometric function can be seen in [8]. In Ref.[8], it is also shown that 
the confluent Heun equation can be normalized to constitute a group of orthogonal com- 
plete functions and the confluent Heun equation also admits quasipolynomial solutions for 
particular values of the parameters. 
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Since D = in our case, it follows from the three-term recurrence relation that 
H{D, A, B,C,h] z) is a polynomial solution if Wi{h) = 0, where Wi stands for a poly- 
nomial of degree 1 in h. Namely , there is only one eigenvalue hi for h such that Wi{hi) — 
(i.e. A = 0). 

V. REDUCTION OF THE ANGULAR EQUATION INTO A SET OF LINEAR 
FIRST ORDER DIFFERENTIAL EQUATIONS: THE CASE WITH MASS 

To complete our analysis of the angular equation, we need to discuss the angular equation 
for the Dirac particles with mass. 

The angular equations (20,21) can be rewritten in the following forms 



LA3 + ifipVl + cos2 9Ai = iXAi (48) 

L+Ai + iUpVl + cos2 eA^ = iXAs (49) 

With substitutions 

A^{e) = {Ao{e) + iBo (9)) ef (50) 

A3{9) = {Mo{9) + iNo (9)) ef ^1^'^^ (51) 
we can transform equations (48,49) into a set of first order differential equations 
M^(^) + GMo{9) - ^iVo (0) = - (a + i^pVl + cos^9) (9) 

N!,{9) + GNo{9) + {9) = (A - /x^Vl + cos^ Aq {9) 

A'o{9) - GAo{9) - (9) = - (A + //^Vl + cos^^) No (9) (52) 

B'^{9) - GBo{9) + {9} = (a - /x^Vl + cos^^) Mq (9) 
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Introducing a new variable x — cos 9 and with the further substitutions 



Mo{e)^-{mo{e) + aQ{e)) 



No{e)^-{no{9) + boie)) 



MO) = o K(^) - ao{9)) 



Bo{e) = - {no{e) - bo{e)) 



(53) 



we may obtain the final form of the set as linear first order differential equations 



mo(x) + aiao(x) + (0:2 — CI3) no{x) = 



ao(a;) + aimo{x) + {a4 + a^) bo{x) = 



rio(a;) + aibo{x) - (q!2 + eta) rno{x) = 



(54) 



6q(x) + aino{x) — (0:4 — 0:3) ao{x) ~ 



where 



ai — — 



m (1 + x^) 



1 



q;2 



A 



2(1 + ^2) 713- 



as = 



Q;2 = 



A 



(55) 



2 (1 + x2) vr^r^ 

Although system (54) does not seem to be solved analytically, one may develop an ap- 
propriate numerical technique to study it. In the literature, there may exist such interesting 
systems which are more or less of this type. 



VI. THE REDUCTION OF DIRAC EQUATION TO ONE-DIMENSIONAL 
SCHRODINGER-TYPE EQUATION WITH A CONSERVED CURRENT 



It is possible to get more compact forms the axial equations (18,19) by introducing the 
scalings 
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g,{y) = Z,{y){l + y'y~' (56) 

^2(2/) = ^2(2/)(l + y')"^ (57) 

and applying the coordinate transformation y — iaxiu , the axial equations take the 
forms 

Z[ (u) - iaZi (u) = -iXXZ2 (u) (58) 

Z!^{u) + iaZ2{u) = iXXZi{u) (59) 



where X — —ri= 



Letting 



cos 



Z,(u) = '-P'f") - -P^f") (60) 



Z,(u) = '^'M + ^'W (61) 
we can combine equations (58,59) to give 

P[{u) = -E+P2{u) (62) 

P^iu) = E_Pi{u) (63) 

where 

E+^a + XX (64) 

E_^a-XX (65) 

Decoupling is attained, by introducing 

Pi{u) = ^/E^T{u) (66) 
16 



P2{u) = ^/EIs{u) 
where we obtain a pair of one-dimensional Schrodinger-type equations 

T"(u) + ViT(u) = 



(67) 



(68) 



S"{u) + V2S{u) = 



(69) 



with the potentials 



) + 



AX 



(a (l + 2y') + XX') 



(70) 



2El 




2El 



XX 



{a (l + 2y') - XX') 



(71) 



One can easily observe that for y ±oo the potentials diverge. This result stems from 
the fact that our space-time is not asymptotically flat. 

To examine the existence of the superradiance, one may consider the conserved net 
current of Dirac particles [1]. In other words, the rate . at which particles falling 
through the horizon per unit time, which must be negative for the superradiance to occur 



where g is the determinant of the spacetime metric and is the axial component of the 
neutrino particle current. We recall from metric (7) that we have 



It is clear from transformations (3,6) that the horizon of metric (7) corresponds to y 
(— oo). In other words, integral (72) is taken over y (— oo). 
In the more standard spinor formalism, is introduced as [1] 




(72) 




(73) 




(74) 
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where 



'AB' 



1 

72 



V 



-V 



(75) 



In this notation, the basic spinors defined by P and Q correspond to [6], 



(76) 



We evaluate as 



'1+1/2 



(77) 



Assuming that the angular functions Ai {6) and {9) are normalized to unity, the 
integral in equation (72) yields that 

J ^/^J^dedip = 3.2467r ([Zsl^ - \Zif) (78) 

Prom equations (60), (61), (62) and (63), we successively find 

(P1P2 - P2P1) (79) 



\Z2\^-\Zif 



2 



2E. 



Pi, Pi 



where 



Pi, Pi 



is the Wronskian. 



Therefore, in order to check the existence of the superradiance, it will suffice to seek the 
solution for Pi at the horizon. 

The reality that the potentials Vi and V2 become infinite both at the horizon and y — > 
00, leads us to think the problem as a problem of particle in the infinite potential well. 
Since the particles are bound inside the well, the principal physical fact requires that the 
solutions of the wave equations (68,69) must be identically zero at the walls (the horizon 
and ?/ — > 00). Clearly, the Wronskian vanishes at the horizon and it follows that the number 
of particles exiting the horizon per unit time is zero. Consequently, similar to the general 
Kerr background [1,3], there is also no superradiance in the extreme Kerr throat geometry. 
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VII. CONCLUSION 



Our aim in this paper was to do more than separating the Dirac equation in a sector 
of Kerr, namely the extremal Kerr throat geometry and obtain exact solutions if possible. 
This premise has mostly been accomplished and they definitely will contribute to the wave 
mechanical aspects of spin-| particles prior inf ailing into the extreme Kerr BH. 

In the general Kerr background the radial Dirac equation was the harder part to be 
tackled compared with the angular part [1]. In the present problem of the extremal Kerr 
throat we have the opposite case: the axial part poses no more difficulty than the angular 
part does. For the massless case, we overcome the difficulty and attain exact solution in 
terms of Heun polynomials. Inclusion of mass prevents this reduction and as a result we 
are unable to express the angular equation in terms of a set of known equations. This part 
of the problem can be handled numerically. Alternatively, the angular equation casts into 
a pair of Schrodinger-type equations. Unlike the scalar field case Dirac fields exhibit no 
superradiance. The charge coupling of a Dirac particle to an extremal Kerr-Newman BH in 
its near horizon limit may reveal more information compared to the present case. This is 
the next stage of study that interests us. 
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